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ABSTRACT 

Context. Dark cloud chemical models usually predict large amounts of O2, often above observational limits. 

Aims. We investigate the reason for this discrepancy from a theoretical point of view, inspired by the studies of Jenkins and Whittet 
on oxygen depletion. 

Methods. We use the gas-grain code Nautilus with an up-to-date gas-phase network to study the sensitivity of the molecular oxygen 
abundance to the oxygen elemental abundance. We use the rate coefficient for the reaction O + OH at 10 K recommended by the 
KIDA (Kinetic Database for Astrochemistry) experts. 

Results. The updates of rate coefficients and branching ratios of the reactions of our gas-phase chemical network, especially N + CN 
and Hj + O, have changed the model sensitivity to the oxygen elemental abundance. In addition, the gas-phase abundances calculated 
with our gas-grain model are less sensitive to the elemental C/0 ratio than those computed with a pure gas-phase model. The grain 
surface chemistry plays the role of a buffer absorbing most of the extra carbon. Finally, to reproduce the low abundance of molecular 
oxygen observed in dark clouds at all times, we need an oxygen elemental abundance smaller than 1.6 x lO"''. 
Conclusions. The chemistry of molecular oxygen in dense clouds is quite sensitive to model parameters that are not necessarily well 
constrained. That O2 abundance may be sensitive to nitrogen chemistry is an indication of the complexity of interstellar chemistry. 

Key words. Astrochemistry - ISM: abundances - ISM: molecules - ISM: individual objects: L134N - ISM: individual objects: 
TMC-1 
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1. Introduction Many explanations have been proposed to reconcile obser- 

vations and models. From a chemical modeling point of view. 

Chemical models predict that gas-phase oxygen should be pure gas-phase chemical models can ex plain the observed uppe r 

mainly in the form of O2 an d CO in t he cold inte rstellar limits for cloud s younger than 10^ yr (jWakelam et alj|2006al). 

medium (see for instance Wak elam et aP |2006a; Oua n et al.l IViti et al.| (|2001|) explore d the possibility that c hemical models 

5008). Since the 1980's, there have been searches for O2 in the can display bistabilities (|Le Bourlot eialj|l993l). In the param- 

inte rstellar mediu m using ground-based and space telescopes eter space where bistability exists, one of the solutions is char- 

(see iPagani et all 12003. and references therein). First, analyses acterized by a very low abundance of O2. In this solution how- 

of data from the SWAS satellite gave an upper limit of about ever, all molecular s pecie abundances are v ery small, including 

10 in dense clouds (Goldsmith et al. 2000). The ODIN satel- the CO abundance (IWakelam et al.ll2006bh. which is not what 

lite also initially gave negative results (Pagani et al.ll2003i) with is observed. As another possibility, the effect of uncertainties in 

improved upper limits o f ~ (1 - 2) x 10'^ for nine sources, the rate coefficient of the main reactio n of product ion of O2 (O 

However, a reanalysis bv lLarsson etaH (l2007l) . using more pre- + OH ^ O2 + H) was explored by |Quan et al.| (12008). One 

cise knowledge of the telescope behavior, resulted in a detection needs however to decrease this rate coefficient by a considerable 

of O2 in p Ophiuchi cloud, with a beam-averaged abundance amount to modify the predicted abundance of O2. 

of 5 X 10"^ relative to H2. Using ground-b ased observations of In the presence of dust, molecular oxygen in the gas phase 

'^O'-^O and C'**0 lines. iLiseau et all (12010 ) argue that the emit- can be adsorbed onto grain surfaces. Adsorption of O2 onto 

ting region may be much smaller than the beam of ODIN thus dust grains is insufficient in itself to lower the O2 abundance 

the O2 abundance could be larger by one or two orders of magni- after 10^ yr, because of the balance with thermal evapora- 

tude. Regardless of the exact numbers, the sparsity of O2 detec- tion and cosmic -ray-induced desorption. However, the adsorbed 

tions in the various target molecular clouds is an indication that O2 molecule can be successively hydrogenated to form HO2 

this molecule may not be a reservoir of oxygen. and H2O2. Then H2O2 reacts with H to form water (see also 
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iRoberts & Herbsjl2002h . which is more difficult to release from 
the grains. This surface chemistry allows one to decrease the 
abundance of O2 after 10^ years, but a peak in the abun- 
dance larger than observational limits remains between 10^ and 
10^ years. Unfortunately, this time range encompasses the "age" 
of cold cores determined by the comparison between large net- 
work chemical models and observations of more than 30 species 
in TMC -1 (CP) and L134N (N) ( see Wakelam et al. 2006a; 
iGarrod et al. 2007; Smi th et al.l2004 . lBergin et al.l(f2000l) found 
a small abundance of O2 in the gas at all times, in agreement with 
observations, using a simplified version of grain surface chem- 
istry assuming conversion of O to H2O and C to CH4, and choos- 
ing the branching ratio of the reaction HaO^ + e" — > H + H2O 
that equals 0.33 to reproduce the H2O abundance. Experimental 
measurements from lJensen e t al.'feOOO) showed that this branch- 
ing ratio is 0.25. Finally, flollenbach et al. (2009) studied the 
chemistry of O2 as a function of the depth in molecular clouds 
using a one-dimensional steady-state PDR model. The authors 
studied the influence of many parameters and found that the 
O2 abundance peaks at Av between 4 and 6 and that molecules 
would be strongly depleted at larger Av. 

In this paper, we revisit the question of the O2 abundance in 
dark clouds using the gas-grain model Nautilus with the most 
recent gas-phase network from the KIDA database and new in- 
sight i nto the oxygen elemental abundances provided by JenkinsI 
(l2009l) . In the next two sections, we introduce the problem of the 
choice of oxygen elemental abundance for dense cloud chemical 
modeling and describe our chemical model. The results of our 
simulations and comparisons with observations in the two dark 
clouds L134N and TMC-1 are shown in section |4] We present 
our conclusions about our work in the last section. 

2. Carbon and oxygen depletion 

The depletion of oxygen related to cosmic refe rence values 
has be en investigated for many years. For instance. iMever et al.l 
(Il998h concluded that there was no evidence for density- 
dependent oxygen depletion from the gas phase based on the 
analysis of their observations of 13 stars using the Goddard 
High Resolution Spectrograph (GHRS) onboard the Hubble 
Space Telescope (HST). They thus proposed a mean interstel- 
lar gas-ph ase oxygen abundance (O/H) of 3.19 x 10"^. A few 
years later. Cartle dge et aP fc004) found a weak correlation be- 
tween the oxygen gas-phase abundance and the density of the 
clouds, analyzing a larger sample of sources observed with the 
Space Telescope Imaging Spectrograph (STIS) of HST and Far 
Ultraviolet Spectroscopic Explorer (FUSE ). At densities below 
1 cm"^, they found an O/H abundance of 3.9 x 10"^ and an abun- 
dance of 2.8 X lO""* at densities above 1 cm"-'. 

A comprehensive study of the depletion problem bv lJenkiiisI 
(l2009 h brought additional findings to light. Jenkins re-analyzed 
archival data of atomic lines in more than two hundred lines of 
sight. The main result of his analysis is that all elements (ex- 
cept nitrogen), and even sulphur, seem to deplete with density. 
The depletion factor for carbon is much smaller than for oxy- 
gen (but mostly less certain) so that one expects an increase in 
the C/O gas-phase elemental ratio with the density of the cloud. 
The main mechanism capable of explaining this depletion is ac- 
cretion onto grains. The depletion of oxygen, however, cannot 
be accounted for by the formation of silicates and oxides since 
the oxygen depletion appears to be larger than that of Mg + Si 
-I- Fe. In the diffuse medium, densities are so low that the col- 
lision probability between grains and gas-phase species cannot 
explain the depletion. It is usually assumed that depletion occurs 



in denser regions and that the mass exchange bet ween the dens e 
and diffuse medium explains these observations (lDrainelll990l) . 
In a sense, the observation of the elemental depletion as a func- 
tion of the cloud density traces the survival of refractory com- 
pounds to the UV radiation in the ISM. This may indicate that 
depletion of the elements is stronger in denser regions. There is 
indirect evidence of this high depletion, such as the small abun- 
dance of SiO in dense clouds (ZiujQ/s et al. 1989) compared to 
that of cosmic silicon. To reproduce observatio ns of gas-phase 
molecules in dense clouds. lOraedel et aH (Il982h proposed to use 
a "depletion factor" of ten on observed atomic abundances in 
diffuse clouds for all elements except He, C, N, and O. These 
elemental abundances constitute the largely used "low metal" el- 
emental abundances. The value of the "depletion factor" is how- 
ever poorly constrained and remains a free parameter in chem- 
ical models. A sensitivity analysis by Wakelam etal. (201()3) 
show that this is one of the most important parameters. 

Carbon and oxygen are usually assumed not to display addi- 
tional depletion compared to diffuse clouds but this view is now 
changing. Additional depletion of oxygen in an organic refrac- 
tory component of interstellar dust seems to be the most proba- 
ble way o f accounting f or all oxygen in dense environments as 
shown bv lWhittetl (12010) . Following fresh insight from Jenkins 
about oxygen depletion, IWh ittet compiled an inventory of the 
different forms of O-bearing species as a function of cloud den- 
sity from an observational point of view (see his Fig. 3). In dif- 
fuse clouds, oxygen is partly depleted in silicates and oxides 
(this is the fraction of cosmic oxygen already depleted in the dif- 
fuse gas), in atomic form (which is observed), and in an unidenti- 
fied form (i.e. "unidentified depleted oxygen" called UDO, pos- 
sibly organic refractory compounds). As the density increases, 
the fraction of UDO increases whereas that of atomic oxygen 
decreases. At densities high enough for chemistry to be effi- 
cient (atomic oxygen then cannot be observed anymore), oxy- 
gen goes into gas-phase CO and ice compounds (mainly CO, 
CO2 and H2O). At hydrogen densities of around 1000 cm"^, 
28% of oxygen would be in gas-phase CO and icy species, 19% 
would be in refractory silicate and oxide forms, and presumably 
49% would in UDO. The net result of this is that only 32% 
of the oxygen would be available for the chemistry. We note 
that the missing oxygen in the d ust phase has been debated in 
IVoshchinnikov & Henningl(l2010l) . 

3. The model 

3.1. Nautilus 

We used the Nautilus chemical model described in lHersant et al.l 
t2009) . The model solves the kinetic equations for gas-phase and 
grain surface chemistries. Details of the physical and chemical 
processes included in the model are given in a benchmark pa- 
per by Semenov et al. (2010). We use typical dense cloud con- 
ditions: a gas and dust temperature of 10 K, an H density of 
2x10^ cm"^, a visual extinction of 10, and a cosmic -ray ioniza- 
tion rate of 1.3x10 s '. A single grain size of 0. 1 is u sed to 
compute adsorption rates following lHasegawa et al.l(ll992 ). The 
cross-section per H nucleus used is 5.72 X lO"^-^ cm^ and we as- 
sumed thermal velocities. The adsorption energy used for O2 is 
1000 K (Herma Cuppen private communication). As cosmic rays 
can penetrate deep into grains, we do not restrict the evaporation 
of species by cosmic rays to the surface layer of molecules. The 
species are assumed to be initially in an atomic form as in dif- 
fuse clouds except for hydrogen, which is converted entirely into 
H2. Elements with an ionization potential below the maximum 
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energy of ambient UV photons (13.6 eV, the ionization energy 
of H atoms) are initially in a singly ionized state, i.e., C, S, Si, 
Fe, Na, Mg, CI, and P. 

3.2. Chemical network 



Table 1. Elemental abundances (/H) 



The chemical network, adapted from iGarrod et al.l (l2007l) . in- 
cludes 6142 reactions, of which 4394 are pure gas-phase re- 
actions and 1748 are grain-surface and gas-grain interaction 
reactions. The model follows the chemistry of 458 gas-phase 
species and 196 species on grains. The gas-phase network 
has been updated according to the recommendations from the 
experts of the KIDA databas^l KIDA, for Kinetic Database 
for Astrochemistry, is a recently opened online database of 
gas-phase reactions of interest for astrochemical (interstellar 
medium and planetary atmospheres) studies (see Wakelam 2010; 
IWakelam et al.l l2010bl) . Recommendations on rate coefficients 
by experts in physico-chemistry are given for key reactions. 

The rate coefficient of the reaction O 4 - OH — > O2 + H is 
a subject of debate (see O uan et al.ll2008l) . The KIDA experts 
suggest that the rate coefficient at 10 K is between 2 x 10"" and 
8 X 10~" cm^ molecule"' s"' (see section [3.2. 11 1. In the rest of 
the paper, we use the lower limit for this rate coefficient. Results 
using the upper limit are discussed in section 1431 



3.2.1 . Experimental and theoretical determination of tine rate 
coefficient of tlie reaction O + OH — > O2 + H 



reaction has attracted considerable 



experimental a t tentio n (iHow ard & Sr mthl 1980i 
1980!; iBrune et al. 1983!; 'Smith & 



1981 



The study of this 
menti 

Lewis & Watsoii' ^1980!; 'Brune et al. 1983!; 'Smith & Stewarl 
1994t [Robertson & Smith 2002, 2006; Carty et al. 2006), 
and there have also been a large numb er of theo retical 
studies using a variet y of meth o ds ( , Harding et al.l 
iTroe & Ushakovl 120011: IXu et al.1 l2007t iLinetail 



Liaueetal.1 120091: lO uemener et al.' '2009t iMaer^oiz et al.1 
2004t lJorfietal.l 1200 9; Li et al. 2010). The experimental 



rate constant is well determined between 140 and 300 K 
decreasing from 7 x 10"" cm-' molecule"' s"' at 140 K to 
3 X 10"" cm^ molecule"' s"' at 300 K. Between 40 and 
140 K, the reaction has been studied in a CRESU (Cinetique 
de Reaction en Ecoulement Supersonique Uniforme) ap- 
paratus (ICartv et al.l 120061) leading to a value of around 
3.5(+1.0) X 10"" cm^ molecule"' s"' which however has 
large uncertainties. Quasi-classical trajectory calculations 
give good agreement with experiment between 300 K and 
3000 K CT roe & Ushakov 2001) and between 40 K and 140 K 
dJorfi et al.l i2008h but the relatively good agreement at low 
temperature may be fortuitous. The reaction, which proceeds 
through a relatively long-lived HO2 complex, should be 
amenable to a statistical treatment and the statistical adiabatic 
channel model shoul d be appropriate (Harding et al. 2000; 
iTroe & Ushakovll2001h but has to deal with dynamical barriers. 
Surface hopping forward and backward between adiabatic 
channel potentials on several electronic states allows the system 
to avoid the dynamical bottleneck to some extent leading to 
a marked increase in the rate consta nt around 50 K. Th e rate 
constant calculated by this method ( iHarding et al.l I2OO0I) may 
be considered as an upper value. Time-dependent wave packet 
methods are unsuitable for the low temperature regime leading 
to an unreliably low rate that is constant at low temperature 
(IXu et al.l I2OO7L iLin et al.1 [2008.) . Time-independent quantum 



Element 


Abundance 


Element 


Abundance 


He 


0.09" 


N 




C 


1.7(-4)'' 





(3.3-2.8-1.8-1.4)(-4)^ 


s 


8(-8)'' 


Si 


8(-9)'' 


Fe 


3(-9)'' 


Na 


2(-9Y 


Mg 


7(-9)'' 


CI 


H-9Y 


P 


2(-10)'' 







Notes. " see IWakelam & Herbsll dlOOg) * from I JenkinsI j2009l) . ' see 



text Low Metal elemental abundances 



mechanical calculations, supposed to be the more accurate at 
low temperature, have been applied to this system leading to 
value around 4 x 10"" cm"* molecule"' s"' at 10 K (Liaue et ajj 
2009; Ouem ener et al.l|2009) . However, they neglect spin orbit 
coupling and electronic fine structure of O and OH, as well 
as surface hopping dynamics between the ground and excited 
potentials at large O-OH distances. Among these effects, surface 
hopping dynamics is suspected to be important at low tempera- 
ture CMaergoiz et al. 200 4; Li et al. 2010) and the rate constant 
may be as high as 8 x 10"" cm^ molecule"' s"' at 10 K. Taking 
into account the various calculations and measurements, the 
KIDA experts recommend that the rate coefficient at 10 K is 
between 2 x 10"" and 8 x 10"" cm^ molecule"' s"'. 



3.3. Elemental abundances 

For each element, we define the elemental abundance as the ratio 
of the number of nuclei both in the gas and on the dust grains to 
the total number of H nuclei. This excludes the nuclei locked in 
the refractory part of the grains. With this definition, the elemen- 
tal abundance would be the "gas phase elemental abundance" at 
high temperature, when all ices are sublimated. Since the amount 
of the oxygen depletion is only based on indirect measurements, 
we consider four values for the oxygen elemental abundance : 
1) 3.3 X 10"'': a "low depletion" case, using the gas-phase abun- 
dance observed in t he diffuse cloud ^Oph, frequently used as a 
reference derived by 'Jenkins (|2009|) (t his is also the mean abun- 
dance observed by Meyer et al.l (ll998l) ): 2) 2.4 x lO""*: an "inter- 



http://kida.obs. u-bordeaux 1 .fr 



mediate depletion"; 3) 1.8 x 10""*, following lWhittell (i20 1 Ot) : this 
value is very close to the "low metal" case; 4) 1.4 x I O"'^: a "high 
depletion" case determined from Fig. 16 of jjenkinsl (l2009l) with 
an extrapolation to a density of 2 x 10^ cm"^. 

For simplicity, we use the "low metal" elemental abundances 
for the other elements except for He, C, and N. Carbon is in a 
similar situation as oxygen in the sense that no additional deple- 
tion is usually assumed, although Jenkins (2009) showed that it 
does exist the tre nd being much weaker and less robust. From 
the observations, iJenkinsI (12009.) determined an atomic carbon 
abundance of 1.7 x 10"'* in ^Oph. We use this value, as the ex- 
trapolation to a density of 2 x lO'' cm"^ using Jenkins' relations 
would only reduce the carbon abundance to 1 .2 x 10"'*. The C/O 
elemental ratio corresponding to our four models are 0.5, 0.7, 
0.9, and 1.2. N itrogen is the only element that was not found by 
iJenkinsI (|2009|) to deplete with density. Although it may be an 
observational bias as argued by Jenkins, we use the gas-phase 
abundance observed in ^Oph of 6 .2 x 10"^. The helium ab un- 
dance is assumed to be 0.09 (see IWakelam & Herbstlllool for 
discussion). All elemental abundances are listed in Table[T] 
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Fig. 1. O2 abundance (/H) as a function of time for four values 
of oxygen depletion (top) and as a function of oxygen elemental 
abundance for four ages (bottom). 



4. Results of the chemical simulations 

4.1. Computed abundances 

Figure [T] shows the computed abundances of O2 in the gas phase 
as a function of time for the four elemental abundances of oxy- 
gen described in section |33] (left panel) and as a function of the 
oxygen elemental abundance for different times (right panel). A 
decrease in the oxygen elemental abundance produces a general 
decrease in the O2 abundances at any time, although the effect 
is stronger between 3x10^ and 2x10^ yr Assuming that O2 
has been searched for in a variety of clouds with ages across this 
range, the non-detection of O2 with abundances above 10"^ com- 
pared to total hydrogen would require the elemental abundance 
of oxygen to be smaller than 1.6 x 10""^. 



Figure |2] shows the computed abundances of a selection of 
gas phase species as a function of time, for the extreme "low" 
and "high depletion" cases. As expected, the abundances of 
carbon-rich species (such as cyanopolyynes) are higher in the 
"high oxygen depletion" case (higher C/O ratio), whereas the 
abundance of O-bearing species is lower The various O-bearing 
species are however not influenced to the same extent. CO, OH, 
and H2O are changed only slightly. O2 and SO2 are lower by 
two orders of magnitude and more than one order of magnitude, 
respectively, at the peak abundance (~ 4 x 10^ yr). 

In our four cases, the elemental C/O ratio varies over a large 
range. Pure gas-phase chemical models w ould be very sensitive 
to these variations: for example, Wakelam et al.l (l2010a ) show 
that the HC7N abundance can be modified by four orders of mag- 
nitude at 10^ yr when C/O goes from 0.7 to 1. In our present 
study and at 10^ yr, the HC7N abundance does not depend much 
on the different values of C/O that we adopted, as can be seen 
in Fig. |2] When the C/O ratio increases, this occurs because the 
available C is mainly u sed to form C-rich m olecules, e.g. C„H„,, 
on the grains (see also iGarrod et al.|[2007l) . However, the mod- 
eled abundances are not only sensitive to the elemental C/O ratio 
but also to the elemental abundances themselves. Increasing both 
elemental abundances by a factor of two would for instance in- 
crease the CO gas-phase abundance by two orders of magnitude 
at 5 X 10* yr in all the models, but the O2 abundance at the peak 
(~ 4 X 10^ yr) remains unchanged. In this dense cloud modeling, 
the CO abundance in the gas phase decreases after a few 10^ yr 
because CO sticks onto grains and is then hydrogenated to form 
H2CO and CH3OH. When the C and O elemental abundances are 
increased, the gas phase CO abundance, and as a consequence 
the solid CO abundance, increases and takes longer to decrease 
its abundance. If we allow the system to evolve up to 10^ yr, the 
CO gas phase abundance is only two times larger than the one 
obtained with our previous elemental abundances. Decreasing C 
and O elemental abundance by a factor of two would slightly 
decrease the O2 abundance. 

We started our chemistry assuming that all species were ini- 
tially atomic, except for hydrogen. If we instead assume that all 
carbon is initially in carbon monoxide, the results change drasti- 
cally before 10^ years. Carbon chains are obviously the most af- 
fected species. Using these initial abundances increases the time 
taken to form the carbon-bearing molecules observed in dense 
clouds by a factor of between 10^ and 10* depending on the 
molecule. 

4.2. Effects of the new rate coefficients 

Among the updates of the network, the new values for the 
rate coefficients and branching ratios of the reactions O + C„H 
(»=2,3), O + C?, C^ 4- H2 and particularly O -H H3 (see also 
IWakelam et all 1201 0^ and N -1- CN have changed the model 
sensitivity to the oxygen elemental abundance. In the "high de- 
pletion" case, the new rates result in smaller O2 abundances 
(4 X 10"^ instead of 3.3 x 10"^ at the abundance peak near 
~ 4 X 10^ yr), while this abundance is not affected in the "low 
depletion" case. The O2 abundance is also unexpectedly sensi- 
tive to the reaction N + CN. Decreasing the N + CN rate coef- 
ficient by a factor 3 increases the CN abundance at 3 x 10^ yr 
by a factor of 2.3, and the reaction CN + O2 becomes an effi- 
cient destruction channel reducing the O2 abundance (see also 
discussion about CN + O2 reaction in section l431 l. This shows 
that the chemistry of a relatively simple molecule can be diffi- 
cult to predict, and that unexpected reactions can be important. 
These results underline the importance of using accurate rate co- 
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10"^ 




time (yr) 

Fig. 2. Gas-phase abundances relative to total hydrogen of a se- 
lection of species as a function of time computed for two dif- 
ferent oxygen elemental abundances; 3.3 x lO"'* (grey line) and 
1.4 X 10-4 (black line). 



efficients in all circumstances. Note that using the upper limit in- 
stead of the lower one for O + OH rate coefficient increases the 
abundance peak of O2 by a factor two in the high depletion case, 
and the maximum elemental abundance of oxygen required to 
reproduce O2 abundance in dark cloud is changed to 1 .5 x 10 

4.3. Agreement with observations in dark clouds 

We compared the abundances predicted by our model for dif- 
ferent depletion cases ("high", "intermediate" and "low") with 
observations of gas-phase molecules in two cold clouds TMC-1 
(CP) (the so-called cyanopolyyne peak) and L134N (N) (north 
peak) to check for the effects of oxygen depletion. The observed 
abundances of some relevant species are listed in Table |2] A 
more complete list can be found in iGarrod et al.l (l2007l) . The 



Table 2. Observational values of abundances relative to total hy- 
drogen in TMC-l(CP) and L134N (N) for some species. 



Species 


n(i)/nH" observed 


n(i)/nH'' observed 




(TMC-1 (CP)) 


(L134N (N)) 


O2 


<3.9(-8)'^ 


<8.5(-8)'' 


CO 


4(-5)^ 


4(-5)'- 


OH 


K-iy 


3.8(-8)'- 


H2O 


<3.5(-8)'' 


<1.5(-7)'' 


SO2 


<5(-10)'' 


<8(-10)'' 


CN 


2.5(-9)'- 


4.1(-10)^ 


HCsN 


2{-9y 


5(-ll)^ 


HC7N 


5(-10)f 


l.O(-ll)^ 



Notes. ° aib) = ax 10*; * from'Pagani et al.' (2003*); see'Smi th et al] 
(20 0l); fromlSnell et al.l (2000); ' from Ohishi et aL (1992) ; ^ from 
■Dickens et all ( 120001) ; « from lOhishi & Kaiful ( 119981) 

observational limits on the O2 abundance relative to total hy- 
drog en in TMC-1 and L134N are 3.85 x 10"^ and 8.5 x 10"^ 
(see Pagani et al.lF2003h . respectively, and they are reproduced 
by our models. SO2 is less efficiently produced in the gas-phase, 
which is in closer agreement with the observations, whereas CN 
is now overproduced (CN abundance relative to total hydrogen 
in TMC-1 (CP) and L134N (N) are, resp ectively, 2.5 x lO -*^ 
(see lSmith etaLll20 07) and 4.1 x 10 '° (see lOhishi et"al][T992h ). 
Because of the importance of CN to the destruction of O2 under- 
lined in section |4j2] this excess of CN may artificially destroy 
O2. However, reducing the rate coefficient of CN -1- O2 by a fac- 
tor of 10 (to simulate a lower abundance of CN) does not affect 
the O2 abundance, because other destruction channels take over, 
in particular C -H O2 — > CO + O. Raising the rate coefficient of 
CN + O2 would however continue to lower O2 abundance. 

For a more global view, we define a pseudo-distance D be- 
tween the models and the observations 

D(t) = ^ ^ |log {XJ"' (f)) - log {Xf)\ , (1) 

j 

where X"*' is the observed abundance of species j and X"'"'' (f) 
the abundance of species j computed by the model at the time f, 
and N is the number of observed species, i.e., 53 in TMC-1 and 
42 in L134N. The smaller the value of D, the closer the agree- 
ment. D is displayed for our three main models in Fig. |3] The O2 
abundance predicted by our model for the "low", "intermediate" 
and "high" depletion cases at the time of maximum agreement 
are, respectively, 4 x lO"'', 3 x 10"^ and 8x10"*^ for TMC-1 and 
3 X lQ-\ 7 X 10-^ and 10"'^ for L134N. The elemental O abun- 
dance does not significantly affect the agreement in the case of 
L134N. However, as expected, the agreement in TMC-1 is closer 
in the "high" oxygen depletion case between 2x10^ and 2x10^ 
years. This closer agreement comes from the higher abundance 
of cyanopolyynes during this range of time. 

We must keep in mind that the data obtained by observa- 
tions indicate that the average abundance along the line of sight, 
and the abundances of species varies very considerably from the 
cloud surface to the cloud interior where the photodesorption is 
less important, as pointed out bv .Hollenbach et al.. (.2009) . 

5. Conclusions 

Our study has demonstrated that the low abundance of O2 in dark 
clouds can be explained by gas-grain chemical models based on 
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10^ 10* 10^ 10^ 10' 

time (yr) 

Fig. 3. Parameter D as a function of time for TMC- 1 (black lines) 
and L134N (grey lines) and three different oxygen elemental 
abundances: 3.3 x 10""* (dotted line), 2.4 x 10""* (dashed line), 
and 1.4 X 10"** (solid line). 



the fresh insig ht i nto elemental o xygen depletion obtained by 
iJenkinsI (l2009l) and lWhitte3 (12010 ). High oxygen depletion also 
improves the overall agreement between models and observa- 
tions for other molecules. We have found that gas-grain models 
are less critically sensitive to the C/O ratio than pure gas chem- 
istry. This is fortunate because the elemental C depletion remains 
poorly constrained by current observational studies. This limited 
sensitivity makes our conclusion more robust. 

Besides this important result about the impact of the O de- 
pletion on the O2 abundance, our study has also revealed that 
unexpected reactions may become very significant when the el- 
emental abundances are modified. The example of the N + CN 
reaction is notable in this respect. This may have consequences 
for models using limited reaction networks: while they may be 
valid over a limited range of input parameters (and for a limited 
number of predicted species), they should not be utilized in other 
conditions without re-assessing their performance. 
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